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Abstract 

We consider the quadratic scalar radius of the pion, (r§ and the 
mixed K — tv scalar radius, (rg xw)- With respect to the second, we point out 
that the more recent (post-1974) experimental results in Kis decays imply a 
value, (rg^xn) = 0.31 ± 0.06 fm^, which is about 2 a above estimates based on 
chiral perturbation theory. On the other hand, we show that this value of {r§ 
suggests the existence of a low mass Ktt resonance. With respect to {rs^^^), we 
contest the central value and accuracy of current evaluations, that give (rg = 
0.61 ± 0.04 fm . Based on experiment, we find a robust lower bound of (j's.tt) ^ 
0.70 ± 0.06 fm'^ and a reliable estimate, {r^ n) = 0.75 ± 0.07 fm^, where the error 
bars are attainable. This implies, in particular, that the chiral result for (rg^) 
is 1.4 cr away from experiment. We also comment on implications about the 
chiral parameter I4, very likely substantially larger (and with larger errors) than 
usually assumed. 
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1. Introduction 



The quadratic scalar radius of the pion, (rg and the mixed K — n (quadratic) scalar radius, (rg j^^), are 
quantities of high interest for chiral perturbation theory calculations, or, more generally, for pion physics. 
Using chiral perturbation theory to one loop they can be related to meson masses and decay constants: '^'^1 



7 / \ / 1 + ^-A^ 3a;(l + x) , 



1 

K 



l^|6(2Mi - M-) log M + 9M,2 ^ _ _ (lO + i^) |. 

"(Lib) 

The second can also be expressed in terms of the chiral constant U; to one loop,[^l 

Here fx and M^, are the decay constants and masses of pion and kaon; = 547 MeV is the 

eta particle mass. We take M^r = 139.57 MeV (the charged pion mass), but choose an average kaon mass, 
Mx — 496 MeV. Erom (l.la,b), Gasser and Leutwylerl^l obtain the theoretical predictions 

(^i,if,r>GL = 0.20 ± 0.05 fm^, (1.3a) 

(r§_jGL = 0.55 ±0.15 fm^; (1.3b) 

the errors come from estimated higher order corrections. 

For (rg j^^) we have experimental information from the decays K^^ and K^. For the first the 
world average of the Particle Data Tablest^' is 

Ao = 0.025 ± 0.006 [K^^] (1.4a) 

and (r|^^) = 6X0/ M^. For the four more modern experimental analyses'^' give the numbers^ 

Ao = < 



0.062 ± 0.024 Artemov et al. (1997) 
0.029 ±0.011 Whitman et al. (1980) 



0.019 ± 0.010 Heintze et al. (1977) 
L 0.008 ± 0.097 Braun et al. (1975). 

If we average them, which is permissible since they are compatible within errors, we find 

Ao = 0.027 ± 0.007 [K^], (1.4c) 

a value in perfect agreement with (1.4a), to which it should equal if neglecting isospin breaking effects. We 
compose (1.4a), (1.4c) to get Aq = 0.026 ± 0.005, and find what we will consider the experimental value 
for the form factor: 

(''s,i^^>exp. = 0.312 ± 0.070 fml (1.5) 



^ This is one of the few cases in which the PDT recommend a number difficult to believe. Perhaps influenced by 
very old determinations (pre-1975) they give the average value 0.006 ± 0.007, incompatible both with isospin 
invariance and with the with post-1974 experiments, and which we disregard. 
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On comparing with a dispersive calculation, (1.5) strongly suggests the existence of a low energy S| Ktt 
resonance. On the other hand, the central value in (1.5) lies clearly outside the error bars of the chiral 
theory prediction, (1.3a). 

There is no direct measurement of (rg ^). Donoghue, Gasser ad Leutwyler^^l used the two-channel 
Omnes-Muskhelishvili method and tttt phase shifts to give what is presented as a precise experiment-based 
estimate; Colangelo, Gasser and Leutwyler review it and, with a minor updating, accept it at present: 

(rg,^) = 0.61 ± 0.04 fm^ (1.6) 

It is difficult to believe that the precision and central value in (1.6) hold at the same time. To get these 
numbers, Donoghue et al. use experimental phase shifts for tttt scattering above the KK threshold, where, 
because one does not measure the process KK KK. the set of measurements is incomplete (as proved 
for example in refs. 7) and where, indeed, different fits give totally different eigenphases (necessary to 
perform the Omnes-Muskhelishvili analysis), as may be seen explicitly in ref. 8. Moreover, they neglect 
multipion contributions which, for the electromagnetic form factor of the pion, account for some 6% of the 
full result. 

As a matter of fact, we will give here a new evaluation (which is the main outcome of the present 
note) and will, in particular, present examples of phases which are compatible with experimental informa- 
tion, as well as with all physical requirements at high energy, and for which the corresponding (rg ^) is 
several standard deviations above (1.6). In particular, we find a safe bound, and a reliable estimate: 

(r-g,^) > 0.70 ± 0.06 fm^ (1.7) 

{rU = 0.75 ± 0.07 fm^. (1.8) 

Moreover, we show that the error bars in (1.8) are attainable. 

It thus follows that, also for (rg ^), the errors due to higher orders are underestimated. 

2. The Omnes-Muskhelishvili method for form factors and radii 



We consider the scalar pion form factor, Fs{t), and the mixed scalar form factor fK^it)- We will also 
discuss the electromagnetic form factor of the pion, F„{t). In terms of these, ^ 

Fs{t)^^^Fs{0){l + l{rl,)t}, (2.1) 

and (r^) is the electromagnetic radius of the pion. For F^, current conservation gives i^7r(0) = 1. The 
values of Fs{0), fKni'O) niay be calculated with chiral dynamics, [•^1 but we will not concern ourselves with 
this here. 

Let us denote by F{t) to any of the three form factors in (2.1), and let S{t) be its phase: 

6{t) = aigF{t), t>sth; (2.2) 

Sth is the threshold, 4M^ or {M^^ + Mk)'^, as the case may be. (We note that in (2.2) we do not understand 

the principal value of the argument; the phase has to be taken as varying continuously with t). The Fermi- 
Watson final state interaction theorem implies that, for t < sq (where sq is the energy at which inelastic 

^ We define the form factors by 

{TT(p)\muUu{0) + mddd{0)\TT{p')} = {2tt)-^ Fs{t), 
{w(p)\{ms -mu)us{0)\K{p')) = (27r)-Vxx(t); 
the meson states are normalized to {p\p') = 2poS{p — p'), and t = {p — p')^- 
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channels become nonnegligible), 6{t) equals a corresponding scattering phase. ^ To be precise, 

6{t) — 5i{t) [P wave tttt phase, for F^r]; sq ^ 1-1 GeV^ 
S{t) = 6^\t) [SO wave (S wave with isospin 0) tttt phase, for Fg]; Sq = 4M^ (2.3) 
S{t) = 6Q'^'^\t) [S^ (S wave with isospin i) Kn phase, for /ktv]] sq ~ 1.5^ GeV^. 

We will assume that we know the phases 6i, 5l^\ 5l^^^\ and thus 6{t), for t < sq. 
At large t, the Brodsky-Farrar counting rules'^' imply that 

F(t) — -4^ — -, (2.4) 

^ 't^oo -Hog'^(-i) ^ ^ 

from which it follows that, unless the phase oscillated at infinity, one must have 

d{t) + (2.5) 

t-»+oo [ logt/t J 

In particular, (2.5) implies that 6{oo) = n. For F^, the Jackson-Farrar calculation!^! gives 

F^{t) c l2'KCFflas{-t)/{-t), 

t — ^oo 

hence v — 1; for the other form factors one cannot prove a similar behavior rigorously in QCD, although it 
is likely that v = 1 also here, i is a scale; for the electromagnetic form factor, it is ~ yl^, with A the QCD 
parameter, but its precise value is generally not known. Nevertheless, the feature that the limit 5{oo) has 
to be reached from above, i.e., that at asymptotic energies 5{t) is larger than tt, seems to be general. 

We will use the Omnes-Muskhelishvili method, I^^' with only one channel, to solve for F in terms 
of (5; it will turn out that the two-channel method is neither necessary nor reliable [the last for the reasons 
explained after Eq. (1.6)]. According to it, we have that, given the condition (2.5), the phase determines 
uniquely F: one has 

F{t) = F(0) exp ( - r ds -4^^ I . (2.6) 



s{s - t) 

From this we get a simple sum rule for the square radius (r^) corresponding to F(t): 

/ 2^ 6 /"^ , 5(s) 

{r^) = - ds^. 2.7 

In general, we will split (r^) as follows: 

{r^)^Qj{so) + Q<p{so) + QGiso). (2.8) 
Here Qj is the piece in (2.7) coming from the region where we know S, 

Qj{so) = - r ds^-^. (2.9a) 

is obtained defining an effective phase that interpolates linearly (in t~^) between the values of 6{t) at 
So and oo: we write 

5eff(t)=^+ [5(so)-7r]^, (2.9b) 

and then set 

Q^{s^)^- r ds^-^. (2.9c) 
Finally, Qa corrects for the difference between 5 and des'- 

^ J so S 



^ Actually, the individual relations in (2.3) also hold at any t (even above so) provided inelasticity for the corre- 
sponding paxtial wave is negligible there. 
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Qj, Q<f are known; Qa has to be fitted or estimated. The decomposition (2.8) is equivalent to decomposing 
F as a product. We integrate explicitly 5es and then we can write 



F{t)=F{Qi)J{t)mG{t); 



.\[l-Sl{sa)/^]so/t^ (2.10) 



What we know about G{t) is that G(0) = 1, and that it is analytic except for the cut sq < i < oo. 
The best way to take this into account is by making a conformal mapping of this cut plane into a disk, 
and expand in the conformal variable, z{t): 

= f^-^. (2.na) 

We then write [-'^-'^1 

G{t) = 1+Aq + ciz + C2Z^ + C3Z^ + ■■■ , (2.11b) 

an expansion that will be convergent for all t inside the cut plane. We can implement the condition 
G{0) = 1, order by order, by writing Aq = — [ciZq + C2Zq + c^Zq + ■ ■ ■] , zq = z{t = 0) = —1/3; the 
expansion then reads. 



G{t) = l + ci{z + 1/3) + C2(z2 - 1/9) + C3(z3 + 1/27) + • ■ ■ , (2.12) 

the Ci being free parameters. 

The contributions to the square radius Q^, Qa may be written explicitly in terms of 5{so), Ci as 

So I TT J 350'-^ 

For the electromagnetic form factor of the pion we take, following ref. 11, ,so = 1.1 GeV^. For we can 
fit experimental data and thus find the c^. These data are in fact precise enough to give two terms: I^^l 

ci = 0.38 ± 0.03, C2 = -0.19 ± 0.03 [For F„\. (2.14) 

The squared charge radius is then 

{rl) = 0.435 ± 0.003 fm^. (2.15) 



3. Dispersive evaluation of the square radii 



3.1. The electromagnetic radius of the pion 



We start with a review of the evaluation of {r^), which will serve as a model for the other two. Although 
in ref. 11 the P wave phase was deduced from the experimental values of F^, we here consider it as given, 
for 4M^ < i < 1 GeV^. Taking its value from ref 11 one has (in fm^) 

Qj = 0.195. (3.1a) 

Moreover, 6{so) = Si{so) = 2.70, and hence 

Q<2, = 0.217. (3.1b) 
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So, if we approximated G{t) = 1 (physically, this is approximately equivalent to neglecting inelastic 
channels), we would underestimate the radius, but not by much, as we get 

Qj + Q<p = 0.412, 

which is 6% below the full value of (r^) as given in (2.15). 

It is not easy to guess Qg, although it is easy to imderstand its sign: in (2.5) we have interpolated 
linearly from sq = 1 GeV^, where 5{so) < tt, to S{oo) — tt; that is to say, systematically below the value tt, 
while we know, from (2.5), that S{t) must approach the asymptotic value of n from above. Thus the phase 
5{s) should rise beyond tt, probably around the energy of the resonance p(1450), to reach its asymptotic 
behaviour (above tt) after that. So we expect 



Qg^6 ds > 1.45 GeV. (3.2) 

logs/i 



Indeed, using this formula with reasonable choices of Sas., t, we get a value near the experimentally 

A 9 1 /2 

measured one for Qq; for example, we obtain the exact result, Qg — 0.024, with t ~ 0.3 GeV , Sas. — 
1.8 GeV. 

3.2. The mixed Kn scalar radius 

We first assume the phase S^'^'^^t) to be given, for t^^'^ < 1.5 GeV, by the resonance _ftr*(1430), whose 
properties we take from the PDT.Pl Its mass is = 1412 ± 6MeV, and its width A = 294 ± 23 MeV; 
we neglect its small inelasticity (~ 7%). We write a Breit-Wigner formula for the phase: 

cotd^ '{t) = —(1 - s/MJBo, q = 

and Bo = 2q{M'^)/r^ = 4.15 ± 0.35. We take sq = 1-5^ GeV^ and then we have 

Qj = 0.050 ±0.025, g<p = 0.087 ± 0.001; Qj + Q<z> = 0.137 ± 0.03. (3.3a) 
This means that Qg is large; in fact, on comparing with the experimental value, Eq. (1.5), we find 

Qg = 0.175 ±0.03. (3.3b) 

The corresponding Ci would also be large, Ci = 7.6 

The sum of Qj and Q^ substantially underestimates the value of the mixed scalar square radius: 
the true phase 6{t) of the form factor would have to go on growing a lot before setting to the asymptotic 
regime (2.5). The size of the phase necessary to produce the large Qg required appears excessive. 

An alternate possibility is the existence of a lower energy resonance (or enhancement; we denote it 
by k), beJow the if*(1430), which some analyses suggest,!^^] ^i^h M« ~ 1 GeV and r« = 400 ± 100 MeV. 
In this case, we approximate the low energy phase, s < sq = 1 GeV^, by writing^ 

cot4'/')(i) = ^(1 - s/Ml)B^, B, = 1.8 ±0.5 
zq 

and find 

Qj = 0.070 ± 0.030, Qs, = 0.180 ±0.006; Qj + Q^, = 0.250 ± 0.030, (3.4a) 
which reproduces well the experimental number with a small Qg, compatible with zero: 

Qg ~ 0.06 ± 0.07. (3.4c) 



* The Breit-Wigner parametrization (3.4) for the k should be considered only as an effective one; it is in fact not 
clear that the phase would reach 90° at M„. 
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3.3. The scalcLT radius of the pion: bounding its value 



Next we consider the quadratic scalar radius of the pion. We will, for the SO phase below i^/^ ^ o.96 GeV, 
take the two fits to experimental data in ref. 13: one possibility is 



cot5^"^(s) = 



Ml 



Ml-s 



2k s - \Ml 



Bo =21.04, 



M2 
Bi = 6.62, 



Bo+Bi 



xVd-o.f. =15.7/(19-3); 



,(0) 



= 782 ± 24 MeV; 
= (0.230 ± 0.010) M^^ 



(3.5a) 



Uncorrelated errors are obtained if replacing the Bi by the parameters x, y with 

Bo = y-x; Bi = 6.62 - 2.59.x- : y = 21.04 ± 0.75, x = 0±2.4. 
This will be referred to as 2Bs. Alternatively, we may take 



(3.5b) 



cot 5, 



(0), 



5I/2 Ml 



Ml 



2k s-\M^ 



Ml 



Bo+Bi 



V4m; 



Bo 



xVd-o.f. = 11.1/(19 - 4); = (0.226 ± 0.015) M'^ 

M<, = 806±21, So = 21.91 ± 0.62, Bi = 20.29 ± 1.55, B2 = 22.53 ± 3.48; 

(3.6) 

this we denote by 3Bs. 

Although we think 2Bs to be more close to reality than 3Bs, and although both give very similar 
results, we include 3Bs because it comprises, within its errors, the SO phase shift by Colangelo, Gasser and 
Leutwyler,[^l which these authors present as very precise and incorporating results from chiral dynamics 
(in addition to analyticity and unitarity) . We find, with self-explanatory notation, 

Qj{t^^^ < 0.96 GeV; 2Bs) = 0.452 ± 0.05 fm^ (3.7a) 

Qj{t^/^ < 0.96 GeV; 3Bs) = 0.440 ± 0.05 fm^ (3.7b) 

Between i^/^ = 0.96 GeV and KK threshold, that we take at t^/^ = 0.992 GeV, we use a fit to 
experimental data as given in Eq. (3.8) of ref. 13(a), and get, in fm^, 

Qj(0.96 GeV < < 0.992 GeV) = 0.013 ± 0.002. 

We then find the numbers 

„ _ f 0.465 ± 0.05 [2Bs] 
^•^ ~ \ 0.453 ± 0.05 [3Bs]. 

To calculate we take the value 

(5^°^(4M|.) = 3.14 ±0.52, 

which covers all the experimental determinations, and get — 0.237 ± 0.02. Therefore, we have 
obtained the result 

Qj + Q^ = 0.70 ± 0.06, (3.9) 

and this comprises both cases 2Bs and 3Bs. 

Eq. (3.9) should be interpreted as providing a lower bound on (rg ,^); it assumes that the phase of 
Fs{s) does not increase for s beyond KK threshold, while, as one would deduce from the similar calculation 
of (r^), and as we will see also in the present case, 6{s) should increase somewhat before decreasing to its 
asymptotic value, 6{oo) — n. We have therefore found the result, 

(rg_^) > 0.70 ± 0.06 fm^ (3.10) 



(3.8) 
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FlGURE 1. The 7 = 0, S-wave phase shift corresponding to 2Bs and Eq. (3.12) (continuous hne) and 
the phase of ref. 6 (dotted hne at low energy). Also shown experimental data points included in the 2Bs 
fit: from the Ki4 experiments, and a point from K2n decay; high energy data of Protopopescu et al. (black 
dots), and Hyams et al., Grayer et al. (open circles), and the s-channel solution of Estabrooks and Martin. 
The dotted line at high energy is the line 5 = n; the dashed line is the interpolation corresponding to Ses, 
Eq. (2.9b). The asymptotic phase (5as., with 0.1 GeV^ < t < 0.35 GeV^, is represented by the thick gray line. 

3.4. The scalcir radius of the pion: calculations 

We can get a first estimate of the remaining quantity needed to calculate (rg ^), Qq, by invoking SU(3) 
invariance. If the k is the SU(3) partner of the a, we inded expect ~ 1 GeV. Identifying QQ{Kn) ~ 
QaiT^), and using (3.4), we find an approximate number, 

(r|^) ~ 0.72 ± 0.09 fm^ (3.11) 

A more sophisticated method to get Qg is as follows. As implied by the experimental data on 
TTTT scattering, the inelasticity is compatible with zero (indeed, the central value is almost equal to 
zero) for the SO wave, within experimental errors, in the energy region 1.1 GeV < .si/2 < 1.5 GeV. It thus 
follows that the phase of Fs{s) must be approximately cqu; il to 4°^(s) for 1.1 GeV < t^/^ < i 42 GeV. 

The phases 6^\s), 6{s) will likely not be equal between 0.992 GeV and 1.1 GeV; however, because 
this is a very short range, and the phases are are equal at both endpoints (in the approximation of neglecting 
inelasticity there), it follows that any reasonable interpolation, e.g., a linear interpolation, will give results 
not very different from what one gets by taking, simply, 

6{s) = S^^\s), in the full region, 0.992 GeV < s^/^ < 1.42 GeV. 

The distortion caused by the inelasticity being nonzero just around 1 GeV is negligible, numerically; later 
we will add the estimated error due to above relation being only approximately true. 



-7- 



-F. J. YNDURAIN- 

We take for Sq^\s) the fit to experimental data in ref. 13(a), Eq. (3.8), 

s^/^ > 0.96 GeV; cq = 1.36 ± 0.05, = 0.802 GeV, M/ = 1.32 GeV . 

The corresponding 5q°^ is shown in Fig. 1. We write, choosing the 2Bs fit for the SO wave below KK 
threshold, 

Qj{sq = 1.42^ GeV^) = Qj(0.9922 GeV^) + Qj(0.9922 to 1.42^ GeV^), 

and we have, in fm^, 

gj(0.9922GeV2) =0.465 ±0.05, ^,7(0.992^ to 1.42^ GeV^) = 0.162 ± 0.002; 
Qj{sq = 1.42^ GeV^) = 0.627 ± 0.05. 

We note that the error in Qj(0.99^ to 1.42^^ GeV'^) is only the error coming from Cq in (3.12); the error 
due to neglect of the inelasticity we expect to be much larger, of the order of 10% to 15%. 

In our present approximation, neglecting inelasticity below 1.42 GeV, we have 5(1.42^ GeV'^) = 
4°^(1.422 GeV^) = 5.10 ± 0.03, hence 

g<g(.so = 1.42^ GeV^) = 0.152 ± 0.001. 

Thus, in this calculation, and adding the estimated error due to neglect of inelasticity between 1 GeV and 
1.42 GeV, we find 

<r|^) = 0.78 ± 0.06 (St.) +" 07 (Inelast.) fm^ (3.13) 

compatible with (3.11). Although the central value here is probably displaced upwards (after all, there is 
some inelasticity), so that (3.13) should probably be considered more like an upper bound, we emphasize 
that this value is attainable. Because experimental data are, at rather less than 1 a, compatible with zero 
inelasticity, it follows that any realistic estimate for (rg ^) must have error bars containing the value (3.13). 
This is one of the reasons why a two-channel evaluation is superfluous. 

It is suggestive that, if we take the asymptotic formula (2.5) for &{t), with i between 0.1 GeV'^ and 
0.35 GeV^, then this coincides, on the average and to a 10% accuracy, with the 5Q'\t), 5^s{t) [the second 
as given by (3.12), (2.9b) with sq = 1.42^ GeV^], for between 1.1 GeV and 2 GeV; see Fig. 1. This 
lends additional credence to our calculation (3.13), and it also suggests a different method of evaluation. 
We note that the asymptotic expression, 

'^as.(s) =n{l 



log s/t 



with ^ = 1, intersects the phase given by Eq. (3.12) at s 1.35^ GeV^, for 0.1 GeV^ < f < 0.35 GeV^ We 
can then use (3.12) for s < 1.35^^ GeV^ and the asymptotic expression Jas., with 0.1 GeV'^ <i < 0.35 GeV'^, 
for s > 1.352 GeV^. This gives 

<r|^) = 0.76 ± 0.06 fm^ (3.14) 

i.e., a result almost identical to (3.13). Indeed, any reasonable interpolation between the asymptotic phase 
and the low energy one would give a similar result. 

Taking into account this, as well as the previous results, we get what we consider a reliable value 
for the scalar radius by writing 

(r|^) = 0.75 ± 0.06 fm^ (3.15) 
which encompasses (3.10), (3.11), (3.13) and (3.14). 
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4. Discussion 



We first say a few words about (rg ^^). The central experimental value, (rg j^^)cxp. = 0.31 ± 0.06 fm^, is 
2(7 above the theoretical prediction of Gasser and Leutwyler,Pl 0.20 ± 0.05 fm'^. It would seem that the 
errors were underestimated by a factor ~ 2 by these authors. The experimental number, together with 
our dispersive evaluation, suggest the existence of the k enhancement. ^^^^ 

We next turn to the scalar radius. The experiment-based evaluation of Donoghue, Gasser and 
Leutwyler,[^'^l (rg ^)dgl — 0.61 ± 0.04 fm^, lies below our lower bound, (3.10), and well below our best 

estimate, 0.75 ± 0.07 fm^. To get a value as low as that of these authors, one would have to assume 
that, contrary to indications from tttt scattering, the phase of Fs{t) would decrease above f^^ = 1 GeV 
to about 1/2 of its asymptotic value; and that this continues to hold up to a very high energy.^ This is of 
course a very unlikely behavior and, what is worse for a calculation based upon experimental phase shifts, 
it is incompatible with what one may get within experimental errors, as proved by our calculation (3.13). 
Prom this it follows that the chiral dynamics calculation at one loop,Pl (r§ J = 0.55±0.15 fm^ lies clearly 
below the value suggested by experiment, Eq. (3.15). 

To finish we say a few words on the value of the chiral constant I4, and the connection with tttt 
parameters. We present a few values for I4, all of them, however, using the DGL'^l value for (rg ^): 

Li=4.4±0.3 [BCT]: [4 = 4.2 ± 0.2 [ABT]: 

k = 4.2 ± 1.0 [DFGS]; I4 = ± 0.2 [Best CGL value; with tttt info]. (4.1) 



Here CGL is ref. 6, ABT and BCT are in ref. 15, and DFGS denotes the paper by Descotes et al.I^^l This 
is to be compared to what one gets, at one loop accuracy, from our results here: 

^4= 5.4 ±0.5, with our best estimate, (3.15); 

[4 > 5.1 ± 0.4, with our bound (3.10). ^ " ' 

Two loop corrections to the various form factors have been evaluated in ref. 17; see also ref. 6. If we accepted 
the value given by Colangelo, Gasser and Leutwylerl^l for the higher order correction, bl^ ~ —0.25, we 
would still find that their estimate in (4.1) is too low. This should have implications for the accuracy of 
their description of tttt scattering, as already mentioned by Descotes et al.,[^^l who get, from fits including 
realistic errors for the tttt phase shifts, much larger errors for I4 than the rest. 

In what respects to the connection with low energy tttt scattering parameters, our value here 
for (rg ^) is in reasonable agreement with the D wave scattering lengths deduced in ref. 13, using the 
Proissart-Gribov representation with correct Regge expressions at high energy. Prom the relation 



19 K 
576 ^2^2' 



and using the numbers^ of ref. 13 for oi and the % we find 

(r|^) = 0.83 ± 0.17 fm^ 



It is difficult to point out wiiere lies the failure in the calculation of Donoghue, Gasser and Leutwyler, as it is of 
the "black-box" type. However, a hint is obtained from their statement (p. 356 of ref. 5) that the values of the 
phases above s^''^ = 1.4 GeV do not significantly affect their results. Contrary to this, our explicit calculations 
show that the contributions to (2.7) from energies above 1.4 GeV are large: of 20% for (rg and of 27% for 
the electromagnetic radius, (rj), where one can check the estimate against experiment. They also provide 28% 
of (ri,K.). 

Specifically, we take a^°' = (18.0 ± 0.2) x 10"* ' = (2.2 ± 0.2) x lO""* M'^ and we have improved the 

value of ai by combining the independent determinations of this quantity from F,r and with the Froissart-Gribov 
representation to get ai = (38.3 ± 0.8) x 10"^ M~^. 
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compatible with our best value (3.14). If we had used the relation between (rg ^) and the S wave phase 
shifts, 

we would have obtained somewhat larger numbers, 

(ri = 1.26 ± 0.26 fm^ 1.14 ± 0.21 fm^ 

depending on whether we use the scattering lengths themselves or the value for the combination 2oq^'' — ^Oq^ 
from the dispersive integral in the Olsson sum rule [ref. 13(a), Eq. (4.3)]. These values are a bit more than 
1.5 (J above (3.14): doubtlessly because the higher order corrections are larger for the relation involving 
the S waves scattering lengths. 
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